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Abstract 

The relation between the quark-gluon description of QCD and the hadronic picture is studied up 
to order a s . The analysis of the spin-1 correlators is developed within the large Nc framework. 
Both representations are shown to be equivalent in the euclidean domain, where the Operator 
Product Expansion is valid. By considering different models for the hadronic spectrum at high 
energies, one is able to recover the a s running in the correlators, to fix the p(770) and ai(1260) 
couplings, and to produce a prediction for the values of the condensates. The Operator Prod- 
uct Expansion is improved by the large Nc resonance theory, extending its range of validity. 
Dispersion relations are employed in order to study the minkowskian region and some conve- 
nient sum rules, specially sensitive to the resonance structure of QCD, are worked out. A first 
experimental estimate of these sum rules allows a cross-check of former determinations of the 
QCD parameters and helps to discern and to discard some of the considered hadronical models. 
Finally, the truncated resonance theory and the Minimal Hadronical Approximation arise as a 
natural approach to the full resonance theory, not as a model. 



1 Introduction 



From many evidences, Quantum Chromodynamics (QCD) has been shown to be the proper 
theory for the strong interactions ^HJ- The Operator Product Expansion (OPE) has resulted a 
very powerful and successful instrument to describe the amplitudes in the domain of deep euclidean 
momenta IH El El HI El E3 El ITT] . However, in the low energy region, the theory in terms 
of quarks and gluons becomes highly non perturbative and these degrees of freedom get confined 
within complex hadronic states. Likewise, the extrapolation of the euclidean OPE information to 
the range of minkowskian momenta is highly non-trivial. 

In the large Nc limit -being Nc the number of colours-, QCD suffers large simplifications 
|12l 1131 ITlj . This limit of QCD will be denoted as QCD^ and it turns out to be a very useful tool 
to understand many features in real QCD, providing an alternative power counting to describe the 
hadronic interactions. Taking Nc — ► oo, keeping a s Nc fixed, there exists a systematic expansion 
of the SU(Nc) gauge theory in powers of 1/Nq, which for Nq = 3 provides a good quantitative 
approximation scheme to the hadronic world. Assuming confinement at Nc — > oo, QCDoo is 
equivalent to a theory with an infinite number of hadronic states where the processes are then given 
by the tree-level exchange of an infinite number of resonances. 

In this paper, we study the spin-1 correlators, 

(q^u - q 2 g^) n XY (q 2 ) = i J dx 4 j** < t{j x j y (o)t } ) , (1) 

with X/Y = V, A and the currents = ^uj^u — ^d^d and = ^7^75^ — ^d^j^d. Only 
the sector of light quarks u/d/s will be considered and we will work under the chiral and large Nc 
limits. We will analyse the V+.A and V — A combinations, H LL = Ii vv +Ii AA and R LR = II VV — H AA 
respectively. 

Dispersion relations are nowadays a widely employed method to relate the theoretical OPE 
results in the euclidean domain with the available experimental data ^ImII(£) ea:p in the positive 
energy region |15l 1161 1171 1181 I19j . In Section 2, a pair of alternative sum-rules are presented, 
providing a comparison and cross-check of former dispersive determinations like Laplace or pinched 
sum-rules. 

We introduce first the usual moment integrals A ( ' n \Q 2 ), which give a largest weight to the low 
energy region (t <C Q 2 ), suppressing the high energy range. However, through the introduction 
of the Legendre polynomials into our sum rules, we may build some particular combinations of 
the moment integrals, namely B <k \Q 2 ), which enhance both low and high energies (t <C Q 2 and 
t » Q 2 ) and produce a stronger suppression on the intermediate region. Hence, we are able to use 
at the same time information from the experimental data (low momenta) and perturbative QCD 
(expected to work at t — ► 00). 

On the other hand, we will also consider the average ilmn(z) of the spectral function ^lmn(t) 
through some rational distributions £ , peaked around t ~ z with a given dispersion (A£) ? which 
suppresses the outer regions. This is an analogous procedure to the Gaussian sum-rules ^Hl and, 
in the limit — > 0, the average would recover the value of the amplitude at t = z. The 
advantage of our distributions £ a is that they only depend on the first moment integrals, still under 
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theoretical control; the influence of higher moments is killed. Unfortunately, although one may prove 
that -ImTl(z) follows an OPE-like power behaviour, narrower and narrower distributions require a 
more precise knowledge of the higher dimension condensates and their anomalous dimensions. In 
addition, the appearance of duality violating terms that cannot be analytically expanded around 
Q 2 = —q 2 — ► +00 may yield observable contributions that are dropped off by the OPE |2L)1 I21j . 

In Section 3, former OPE calculations are revisited under the perspective of these sum rules. 
An analysis of the amplitudes in purely perturbative QCD (pQCD) is also performed, with all 
the condensates and duality violating terms set to zero. The range of validity of the OPE and 
pQCD is reduced to the first moments, diverging once we go to higher orders. Matching the 
averaged correlator ilmII M (z) to the experimental one requires accurate information about the 
condensates of high dimension. Nonetheless, in the V + A case, one finds that pQCD seems to work 
fine for energies up to z ~ 1 GeV 2 , pointing out the more reduced impact from the OPE condensates 
in this channel. The phenomenological analysis of the experimental data |22| I23| 1241 125[ I26[ I27j in 
order to determine the OPE parameters (a s and the condensates) is relegated to a next work. An 
alternative derivation seems relevant since there is still some controversy on the values of the higher 
dimension condensates [HI El El EE] • 

In Section 4, we study large Nc QCD and its manifestation into a meson theory with an infinite 
number of narrow- width resonances (RxT <oo) ). The x denotes that our hadron theory must be built 
up chiral invariant in order to ensure the right low energy dynamics |281 1291 1411 149j , although this 
detail is not relevant for the present work. First of all, a resonance theory dual to QCD must recover 
the free-quark logarithm in Il LL (q 2 ) and the 1/Q 2m OPE structure in U LR (q 2 ) [T31 I2H 13U1 1331 1311 1351 
ESI EH EE1 EHj • The novelty of this work is to introduce the conditions required to recover the 0(a s ) 
running in U LL (q 2 ). Reproducing the a 2 logarithmic dependence of the condensate anomalous 
dimensions is, however, a rather complicate problem that goes beyond this work. We will make 
the identification QCD^ = RxT°° since the resonance theory recovers the OPE in the euclidean 
domain providing, in addition, further information of QCD. For instance, Duality violating terms 
exp [— p Q] lacking in the OPE (Q = \J—q 2 ) can be handled and the positive q 2 range becomes 
accessible. 

Once the analysis is taken up to 0(a s ), one is aware that two different energy regimes must be 
considered; the spectral function — Imll(i) will be split into a perturbative part with t greater than 
some separation scale t p ~ 2 GeV 2 , responsible of the pQCD behaviour, and a non-perturbative 
part with t < t p , essential to recover the right OPE 1/Q 2m structure. In the resonance picture, 
a similar splitting is required. The infinite resonance summation in the spectral function is also 
separated into a perturbative and a non-perturbative sub-series. The perturbative sub-series is 
fixed by pQCD, once a model for the asymptotic spectrum of meson masses M 2 is assumed. Due to 
0{a s ) corrections, the parameters of the light resonances in the non-perturbative range may suffer 
important variations with respect to the asymptotic behaviour of the spectrum. They will be fixed 
through a short-distance matching to the OPE. 

In Section 4.3, some available model for the resonance mass spectrum are studied |211 I3U1 
E2 E2 ESI EH E3 EE1 E3, getting a set of predictions for the p(770) and a x (1260) parameters, 
together with the OPE condensates of dimensions four and six in Ii LL {q 2 ) and H LR (q 2 ) respectively. 
Through a five-dimensional model [22], we exemplify how the resonance models implicitly include 
the OPE information together with the duality violating terms. A more exhaustive analysis have 
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been recently done for other models in Ref. |21| . 

In Section 4.4, the Minimal Hadronical Approximation (MHA) at large Nc |4U1 I41j arises 

(oc) 

naturally as a low energy theory of RxT where the infinite series of mesons is truncated. The 
lightest resonance parameters encode the 1/Q 2m information coming from the larger mass states. 
A very successful phenomenology already exists at large N c |HJ EH ESI EH IH3 IHE EH EE] • This 
framework has allowed the developing of robust calculations at next-to-leading order in 1/Nc |49| 
E01EUE1E1IS1E3], achieving a good control of the final state interactions 03 E3 EH1 ES E2B ED 

In Section 5, the dispersion relations developed before are applied to RxT°° ■ Through the 
usual moment integrals ^4 (n) (Q 2 ) we show the equivalence with pQCD and the OPE. The real 
improvement of RxT°° with respect to to the OPE appears manifestly through the B k \Q 2 ) sum 
rules. The physical components £> (fe> (Q 2 ) exp oscillate as k grows, damping off beyond some k, whereas 
the OPE yields a divergent non-oscillating behaviour. The large Nc resonance theory naturally 
reproduces the oscillation although it never vanish since the states own zero-widths. However, 
one finds a pretty good agreement with the phenomenology for the first components, where the 
damping is still not present. This allows considering in Section 6 the averaged amplitudes -ImIL(z) 
and the exploration of the different models for the spectrum. We are actually sensitive to the 
asymptotic behaviour of the mass spectrum M 2 , being some hadronical models more favoured by 
the phenomenology. 

The paper is, therefore, separated in three differentiated parts: In Section 2, we introduce 
the theoretical tools. In Section 3, we revisit the general features of QCD within the OPE and 
pQCD frameworks. Finally, Sections 4, 5 and 6 are devoted to the study of the large Nc resonance 
description and its connection with the experimental data and the OPE. In Section 7, the results 
are summarised and some final conclusions are extracted. 



2 Dispersion relations in QCD correlators 

The perturbative calculation of the vector correlator at lowest order in the a s expansion, 0(a®), 
is provided by the free-quark loop. In dimensional regularization one has: 



n vv (g 2 



N c 
24vr 2 



Aoo + In- 



with the logarithmic ultraviolet divergence AooOu) = <fzrj + 1E- ln4vr, being lE ^ 0.5772 the 
Euler constant and \x the energy scale introduced in the renormalization procedure. 

One useful way to get rid of the renormalization ambiguity and the ultraviolet divergences is 
through the Adler-function [T7| . 



d-A 



+ 0(a s ) 



(2) 



d 

d\nq 



2 U w(<l 2 



(3) 
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with Q 2 = —q 2 . This function carries the whole information of the vector correlator, which can be 
reconstructed through 

/z=~q 2 ^ z 
— A vv (z), (4) 
z =-% z 

once a given renormalization prescription H vv {qo) is provided. 



2.1 Moments A " (Q 2 ) of a correlator 



In general, for any given correlator Il(q 2 ), it is possible to consider a set of more general moment 
integrals with a larger number of derivatives |18j : 



A in) {Q 2 ) = 1 





' d ' 




dq 2 . 



n(<? 2 ) , 



(5) 



where by construction one includes the correlator A 1 (Q 2 ) = Il(q 2 ) and A ( ' (Q 2 ) = A(Q 2 ) the 
usual Adler function. 

These functions are related with the imaginary part of the correlator through the dispersion 
relations 



A (n) (Q 2 ) 



°° dt Q 2n 1 

-Imn(t) 



(6) 



{t + Q 2 ) n+l 7T 

with a large enough number of subtractions so the integral is convergent (n > 1 for the vector and 
axial correlators). 

t-O 2 

Eq. © can be written in a slightly different way through the change of variable x = y^l- 



A (n) (Q 2 ) 



+i dx (i^n 

adj on 
-1 2, 



1. 



-imn 



7T 



Q 



l + x 
1 — X 



(7) 



with n = 1, 2, .. 



The moment integrals of the spin-1 correlators with n > 1 are therefore physical quantities, free 
of ultraviolet divergences; on the contrary to the correlator, the A^ (Q 2 ) are finite and renormal- 
ization scale independent. 



Eq. (J7J) shows that the moment integrals are simply the projections of the function 



a z (x) 



-imn 



7T 



l+x 
1 — X 



(8) 



in the different directions of the non-orthogonal basis of polynomials \p n {x) = i(l — x) n 1 \ of 

I J n=l 

the Hilbert space of real functions L 2 {— 1, +1), with the scalar product ( / | g ) = y <ix /(x) ^(x): 



A (n) (z) = ( p n | a z ) 



(9) 
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2.2 Orthonormal decomposition B W (Q 2 ) of the correlator 



The non-orthogonal basis {pi(x)} is not very convenient in order to recover the absorptive part 
of the correlator. We can rewrite the observable cr z (x) in terms of the orthonormal basis provided 
by the Legendre polynomial Pk{x) (Pi(x) = 1, P2(x) = x...)\ 



9k(x) = (-1) 



k-i 



'2k — 1 



Pk-i(x) = ]T M B k ? Pl (x), 



(10) 



with k = 1, 2, ... and related to the former basis through some given constants M kl . The 

r+l 

vectors of the new basis obey / dx g m {x) g n (x) = & m ,n- This provides for the imaginary part of 
the correlator the spectral decomposition 



k=i 

with the different components given by the projections 



(11) 



B W {z) = (g k \a z ) = J2 M ki A w {z), 

i=i 



(12) 



(fc) 

where the B (z) depend just on the lowest moments. They can be calculated as well through the 



dispersion relation 

r-l 



B (z 



J ^dx g k ( 



x) ■ a z (x) 



dt 



t - z 



t + z 



{z + tf 

The components in the Legendre basis are bounded if the spectral function is finite: 



B (n \z) < yf( a z \a z ) 



Y^\B {k) (z)\ < V2 max{|a 2 (x)|} . 
\ k=i 



-imn(t) . 

7T 

(13) 



(14) 



The difference with other dispersion relations is the use of the Legendre polynomials to pinch the 
dispersive integral. For the moments A (n) (Q 2 ) one employs a weighted distribution which enhances 
the low energy region, decreases at the range of intermediate momenta, and vanishes at t — ► +oo. 

(fc) o 

The distribution for the components B (Q ) is completely different: The dispersion integral is 
enhanced both around t = and t — > +oo, introducing a strong suppression of the integrand at 
intermediate energies around t = Q 2 . Since the experimental data only reaches up to some finite 
energy and local duality is expected to work at very high energies, this procedure allows minimising 
the uncertainties due to the absence of data at intermediate energies. The Legendre polynomials 
allows replacing the lacking data by the pQCD minkowskian amplitude, reducing the impact of 
duality violations in the transition from the experimental data to perturbative QCD. 
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2.3 Spectral function reconstruction 



One extracts several conclusions from Eq. . First to notice is that the expression is an identity 
for any x and z, although a partial knowledge on the moments introduces wrong dependences. 

_____ (fe^ 

The errors in Eq. due to uncertainties on the B (z) are smaller around x = (t = z) 
whereas large fluctuations occur at the extremes of the interval, x = — 1 and x = 1 (t = and 
t = +oo respectively), where the Legendre polynomial reach their absolute maxima and minima, 
-Pfc(±l) = (±l) fc . Thus, the optimal point for Eq. corresponds to x = 0: 

-I oo 

<7*(0) = -Imn(z) = J2 # (fc) W • 9k(0), (15) 

fc=i 



TT 



where one has g k (0) = for k even, and |^(0)| = y/2*=± \P k -i(0)\ = J ^H r( l ( % G 

fcl V 2 + 2 J 

for odd, with the signs provided by (— l) -2- • r(/c) is the Euler Gamma function. 



l 

V2' 



(n) 

The components i3 (z) related to a physical spectral function (which remains finite) become 
smaller and smaller at a certain k since the norm ( a z \ a z ) is bounded, and the series in Eq. (|15|) con- 
verges. For instance, the spectral function corresponding to the vector correlator in the free-quark 
limit (Ii vv {q 2 ) = — ^jfrln-^r) is easily reconstructed from its components B^ v (z) = N £±^? 4,1- 

The eventual knowledge of the components B (z) at all orders in k allows the exact recovering 
of the spectral function at any energy, in particular at q 2 = z. Actually, if one knows a large 

____ 

enough amount of components B (z), such that the remaining terms in the series of Eq. (|15[) 
already converge, then it is possible to give an estimate of -ImII(z). The truncation error would 

(fe) 

be provided by the size of the last components B (z). 

The relation in Eq. (|15|) is stable, i.e., small variations on the B (z) produce tiny fluctuations 
on -ImII(z). This must not be confused with the fact that tiny modifications on the value of the 
correlator at q 2 = —z < may produce (and produces) large instabilities on the tower of components 
and, hence, on the time-like correlator. 



2.4 Averaged amplitudes 

In many situations the description of the amplitude in some energy range may be complicated 
from the theoretical point of view. In these cases, it is sometimes more convenient to consider the 
amplitude averaged through some distribution peaked around a given energy, in the fashion of the 
Gaussian sum rules . 

We will perform the average of the spectral functions ilmll(t) in the :r-space given by the change 
of variable t = z f^~f J > being z the energy of interest. The spectral function is then provided by 
cr z (x) = ^rlmll z ({^f ) ■ The central point x = of the averaging distribution £ a (x) corresponds 
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Figure 1: Distribution £ a (x) for a = 8. Its corresponding components £ ai fc are shown on the right- 
hand-side. The distribution is shown together with the reference data a^ R {x) exp [22] f° r z = 1 GeV 2 . 



to t = z. This mapping of t allows a simpler analysis in terms of the moments. We consider the 
family of distributions 

Uz) = A % , (16) 



' a | 3 

being a > an even number and N a 



r ( 2 a , 



a constant that normalises the distribution 



v^r(|+i)_ 

to 1. This functions are centered at zero (( x = 0) and have dispersion (Ax)? = ^3 . Hence this 
distribution covers the spectral function around t = z within an interval At ~ 2 z Ax. 

Since £ a (x) is a polynomial of degree a, it accepts a decomposition in terms of the orthonormal 
basis of Legendre polynomial {gk(x)}°° ■ 



k=V 

o+l 



£a{x) = ^ k 9k(x) , (17) 



k=l 



where £ ttj i = l/y/2 due to the normalization J dx £ (a?) = 1, and the constant terms £ aj fc with 

even k are zero due to the parity of £ a (x)- This distribution only depends on the first a + 1 
Legendre polynomials. In Fig. (JTJ) we show the distribution and components of £&(x), compared 
to the experimental data a z (x) exp for z = 1 GeV 2 |22| . For the case with a = 8 we have £s,3 = 
2^ t _ 27V2 c _ 2V2E „„ j t _ 7^ 

tr - . &,5 - T43~' »,7 ST and S8.9 - 4862 • 

The mean value of the spectral function is defined through the average 

-ImU(zf = (a z ) u = f 1 dxa z (x)Ux), (18) 
vr J —i 

where the orthonormal decomposition of the correlator in Eq. (|11[) yields 

a+1 

;mIT(z) ° = 

7T 



±ImII(*) €o = iB W (z) . (19) 



k=l 
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Only the first (a + 1) moments are relevant for the amplitude averaged through £ a . This will be 
useful when our control on the high order components is reduced. 

This procedure is analogous to the Gauss- Weirstrass transform of the correlator, where the 
amplitude is averaged through a Gaussian distribution ^5] . One would recover Imn(z) ^lmn(z) 
in the limit when a — ► do. Nonetheless, our theoretical control on the QCD components B (z) gets 
worse as k increases and one needs to go to high enough energies in order to make them reliable. 



3 Perturbative QCD and the operator product expansion 

The Operator Product Expansion in perturbative QCD provides a systematic procedure to com- 
pute the two-point Green functions at any order in a s or operator dimension in the deep euclidean 
regime Q 2 = —q 2 3> h? QCD [HIE! For the spin-1 correlators one has 

U(-Q 2 f PE = «9 (0) ) + E -J^i (20) 

m=2 t « 

where the coefficients ( 0( 2m ) ) are provided by the dimension-(2m) operator in the OPE and they 
depend weakly on the momenta (only through logarithms) [5|. In this work, the term (O (0) ) corre- 
sponds to the identity operator in the OPE and yields the purely perturbative QCD contribution 
(pQCD). The V + A correlator becomes ( C { " ) a -^° U LL (-Q 2 ) free = In ^ in the free quark 

limit whereas for the V—A Green function vanishes for any value of a s (( O ) = n LB (—Q 2 )f ree = 0). 

V-A correlator 

In the case of the V — A correlator (Ji LR = II VV — II AA ) the OPE starts at the dimension six 
operator [S]: 

n L A-Qr FE = (2D 

m=3 

At high euclidean momenta the correlator is driven by the dimension six condensate, 
with (O ) = — / Lira s (qq} 2 at large Nc- We will not consider the anomalous dimensions of the 
condensates ( (2m) ), which will be taken as constants. In this case, one gets for the moments a 
l/Q 2m power structure, 

LR 

„(") s rt 2\° PE ST* °(«.2m) -.i LR (m-l + n)! LR 

Ar(Q) = Wlth a M = (m _iVn' (e W* ( 22 ) 
and a similar thing happens for the components B (Q ), 

k 

with b,. 9 , = V MZ n a, „ . , 23 

(fc,2m) Z_^/ K,n ( n ,2m) ' v ' 

n=l 



OPE 



OO ^ 

(k,2m) 

3 



2m 
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being the constants given by the basis transformation in Eq. (|12|) that relates A^ (Q 2 ) and 

B (Q 2 )- When considering truncated OPE series, both All(Q 2 ) and B ( ^ R (Q 2 ) diverge beyond 
some order n. 

The spectral function, related to the components B^ (Q 2 ) through Eq. (|15|) . formally shows the 
same power behaviour, 

1 LB 

i oo — Imll i oo 

^) = £ ^V^' with ^ ImI C> = ECo • *(°)- ( 24 ) 

m=3 k=l 

Unfortunately, within the OPE the terms 6 (fc 2m) in the coefficients ^Imll } go on growing as k — > oo 
and the summation diverges, preventing the theoretical determination of the spectral function. 

V+A correlator 

Considering the Renormalization Group Equations up to 0(a s ), the V + A correlator 
= H vv + II AA ) shows the structure |Sj 



" 12vr 2 



« 2 2 ft la,(Q 2 ) 



m=2 



with Cf = N 2 C Nc ' anc ^ ^ = — B ~^ n f) provided by the /3-function at lowest 

order, ^jjj^f = P(&s) = ^/?ia s + 0(a 2 ), being a s (n 2 ) the strong running coupling constant. 

For pQCD, with all the condensates set to zero, one finds the moments 

M 2 .vqcd 1 Ng ( 3C F a s (Q 2 ) \ n( n 2frflx\ or\ 
A ll(Q ) = nU^y ^ + 0( - a ^ Q >>' (26) 

which provides the components 

BZ(Q 2 f CD = V2^ 2 (l + ^^1) *M + 0(a 2 s (Q 2 )). (27) 

The perturbative expansion of pQCD in powers of a s produces the spectral function 

1 j rr I \ vQCD \^ K (k) f \ pQCD fn\ Nc f-i , 3C F a s (z)\ 2 
~ImII £i (2) = 2^ B (z) • g k (0) = 1 + — + 0(a fl (z)). (28) 

7T [ 127T \ 4 7T / 

Nevertheless, once the higher dimension operators are taken into account (( ) = -^oi s { G a ^ v G^ u ), 

( Ofgj ) = |vtq s ( ) 2 , ...) the series of B ( ^ L {z) becomes divergent when k —* oo, as it happened before 
in the V — A case. 



3.1 Divergence of the components B w (Q 2 ) within the OPE framework 

Since the moments of the truncated OPE amplitudes diverge at some point, it is important to 
study how and when this divergence occurs. Since we are now interested on dealing with finite 
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spectral functions, the pion pole is removed from the correlators for the analysis in this section: 



n LR {—Q) 



2-. no — 7T 



n M (-Q 2 ) + 



if 



-Q 2 



which yields the moments 



,(«) 



n M (-g 2 



EL 

Q 2 



Q 



(29) 



(30) 



with F ~ = 0.0924 MeV. This allows working with squared integrable spectral functions 



a z (x)' 1 " 77 , owning a convergent series of components £> (fc> (z)^^ 



(fc) 

In order to estimate the physical components B (Q 



components obtained from the experimental ansatz, 
ImII(f) = ImHYt) • 6 (t Q - 1) 



+ 



we show in Figs. © and © the 
to) , (31) 



Imn(g 2 ) pQCD • 0(f 



being ImII(t) e the experimental data with the pion pole removed |22| and to = 3 GeV 2 . In 
order to recover global duality one cannot take an arbitrary matching point |19| I38| , though the 
duality becomes better and better as to is increased and the oscillations in the spectral function 
vanish. Although the experimental data only reach up to t = M 2 in the r-decay experiments j221 
1231 I24j . one knows from the e + e~ experiments that pQCD provides an appropriate description 
of the vector spectral function [2j. Even the r data for the V + A correlator already show this 
regularity for t > 2 GeV 2 . Our experimental ansatz is however less accurate in the V — A channel 
where the fluctuations are wider than in ^ImIl LL (t) and they still remain sizable at the r-threshold. 
Nonetheless, for Q 2 ~ to, the B k (Q 2 ) sum rules suppress the transition region t ~ to as k grows, 
providing this ansatz a first estimate of the physical components. 

We can see in Fig. (j2j) that, although the leading order OPE contribution governs the very 



first components, these eventually diverge. We have taken the value ( O 



(fc) 

order to illustrate the behaviour of the B (Q 



(6) 



-3.8 • 10" 3 GeV 6 in 



jS] . Adding the contribution from some higher 



dimension operators does not solve the problem of the convergence at low energies (Q 2 = 2 GeV 2 ), 
as one can see in the first plot of Fig. Q . We have used the values ( O 



(O 



(10) 



-9.1 • HT 3 GeV 10 , (O 



(8) 



(12) 



6.0 • 10- 3 GeV 8 , 

11.9 • 10" 3 GeV 12 from the review in Ref. 8 , although 
there is still some controversy about the value of the higher dimension condensates [HI El El HOj • 

However, the situation improves drastically once we go to the deep euclidean regime. At 
Q 2 = 10 GeV 2 (second plot in Fig. (©), the addition of higher dimension contributions produces 



appreciable improvements in the convergence of the B\l{Q 2 ) no ^. In Fig. © and next fi gures, all 
the plots are normalized such that -ImII iL (t) 



VQCD t- 



LR 

OO 



1. 



A similar result is obtained for the V + A correlator, where the inclusion of higher 
dimension operators improves slightly the description in the deep euclidean region (Q 2 = 
10 GeV 2 , second plot in Fig. ©). We have considered the chiral limit and used 



the values a s (Q 2 ) = -2vr/ [ft In (Q 2 /A 2 C J 
( ^ ) = lk a s( G^GT ) * 1.2 ■ io- 3 HU, ( O 



With A QCD ~ 



((.) 



L a s (qq] 



0.235 GeV and 
0.85 • 10" 3 0. 



/ 



3 0, 
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Figure 2: Components B LR (Q 2 ) no _ 1T for the V — A correlators provided by the OPE at Q 2 = 2 GeV 2 
and Q 2 = 10 GeV 2 (left and right-hand-side, respectively). It is compared to the components 
derived from the experimental ansatz. From now on, the amplitudes in the plots are normalized 
such that ^ImU LL (t) pQCD * z ^? 1. 




Figure 3: Components B LL {Q 2 ) no _ 7T for the V + A correlators provided by the OPE at Q 2 = 2 GeV 2 
and Q 2 = 10 GeV 2 . 



3.2 Averaged correlators within the OPE 

We have seen that the OPE remains reliable only up to a given finite order k in the moments. 
Through the consideration of a convenient distribution function £ a (a;), we can nevertheless remove 
the influence of the components B (Q ) of order k > a + 1. The averaged OPE amplitude is 
obtained from the calculation of the moments in the euclidean domain and it is compared to the 
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averaged amplitude coming from the experimental ansatz in the minkowskian region. 



The transform of the 1/Q 2m terms, e.g. in the V — A correlator, produces again an OPE- like 
power series: 

LR 1- = LR,U 

(2m) 

Q 2m 7T~ 

,,„ 2 , ^ 

, and where the 



(O"" ) 1 00 -ImIT 

n Lfl (-Q 2 ) = E -tS^ -> b^M* = E " J m) , (32) 



m=3 ^ ' m=|+l 



with the coefficients -Imll, , = (—1)2 (Q ) 

TT (2m) V / > (2m) ' 

condensates (0 LR , ) have been taken as a constant. 

> (2m) ' 



2a +i r(§+|) r(m+f+i) 
v^Fr(a+2)r(§ +i) r(m-§) 



For a given £a(zOj the average of a rational term Tl(—Q ) ~ — ^ s then zero whenever 

2m < a, so the averaged amplitude is never sensitive to the presence of the pion pole. This means 
that, in the short distance region where the V — A correlator is expected to be governed by the 



,LR 



{O _ 
— — term, the averaged spectral function is zero for a > 6. Likewise, ^ImIl LR (z) a accepts the 

same number of Weinberg sum-rules as the original amplitude ilmIl iH (i). Its non-zero experimental 
value for £s( x ) (first plot in Fig. ijlj)) hints that for z < 5 GeV 2 the corrections from (C (2m) ) with 
m > 5 produce rather relevant effects. 

In U LL (— Q 2 ) pQ ° D up to 0(a s ), one finds from Eqs. ()27() and (|2*S)) that for any £ a (x), 

+ 0(a 2 s ). (33) 

It would be interesting to study how this identity is modified at higher orders in a s . The comparison 
of the theoretical V + A averaged correlator to our experimental ansatz is provided on the second 
plot from Fig. (pE|). showing a very good agreement within 0(a 2 ) uncertainties -solid band-. For 



-lmU LL (q 2 ) = -lmU LL (q 2 ) 

7T TT 



the average is independent of (O }, (O ) and (O^), being governed by the pQCD 

distribution. On the contrary to H LR (— Q 2 ), the contributions from the condensates {O ) with 
m > 5 seem to be much more suppressed. 

By considering appropriate distributions £ a (%) ° ne has, in addition, the possibility of recovering 
a s and the vector condensates ( O ) through a direct analysis of the e + e~ experimental mea- 
surements |26l I27| . Nonetheless, this analysis escapes the aim of this article and it is left for a next 
work. The situation in the V — A case seems not so favourable since the energy range of data is 
much more reduced (t < M 2 ~ 3 GeV 2 ), mainly laying within the long distance region where the 
OPE stops being valid. 



4 QCD correlators in resonance theory at large iV^ 

In order to inspect the range of energies where the OPE breaks down it is necessary to add 
some extra ingredients to our theory. One has to add a mechanism that fully reproduces the OPE 
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OPEupto-cO, > and <0, > 

._. OPEupto<O fIfl; > 
OPE up to <0 (12) > 
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pQCD 
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Figure 4: V — A and V + A correlators cr z (x) averaged by the distribution £s(x) 



^ImH(z) ia ). The solid band on the second plot represents the O(aj) uncertainties. 



'«8 



expressions at high euclidean momenta, being able, in addition, to provide the right low energy 
theory (xPT) and the minkowskian description (experimental data, resonance structure...). We 
will see how a theory with explicit resonances suits this picture. 

The exact recovering of perturbative QCD requires the full pile of hadronic states. However, 
the infinite summation of resonance from large Nc QCD is by itself a limit. One of the available 
methods to regularise the series is through an ultraviolet cut-off. Hence, one constructs a set of 
quantum field theories, RxT^ n \ with a finite number n of resonances to describe the observables 
at some given momentum p 11 |351 1361 I37| . QCD^ would be recovered by taking the limit n — > oo 
keeping the external momenta fixed (RxT^[p^} Rjfl<°^ \p»] = QCD 00 [p' 1 }). 



Accommodating RxT^ 00 ^ and QCD at large Nq is not trivial. The first approach must be 
focused on recovering pQCD in the deep euclidean domain. The appearance at Nc — > oo of 
non-perturbative condensates in the OPE produces small corrections into the pQCD correlator at 
Q 2 S> Aq CD - However, due to the non-perturbative effects at large Nc, the smooth pQCD spectral 
function turns into a meromorphic function with an infinite number of real poles. 



It is important to recall the possible existence of duality violating terms. The Green-functions 
( T{ Jn(x) ^(0)} ) are not fully determined by the singularities at x 2 = j^Hi- F° r instance, one may 
have terms of the form gr+s , whose Fourier transform falls off in the deep euclidean as exp (—Qp) 

but becomes oscillating, sm (qp), in the minkowskian region (Q = \/—q 2 and q = y^q 2 )- It is 
actually interesting that some resonance models generate explicitly this kind of duality violating 
terms pointing out the fact that a resonance theory is fully able to produce both the OPE 

and non-OPE QCD singularities. We will derive the OPE condensates in this section by neglecting 
the duality violating terms in the deep euclidean domain although further works should focus on 
the estimate of these terms [2^ . 
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In real world, the resonance poles become complex and get shifted to unphysical Riemann sheets 
due to Dyson- Schwinger resummations at all orders in 1/Nc; once the higher 1/Nc corrections are 
taken into account and the resonance gain their physical non-zero widths the amplitude becomes 
again smooth, as it is found experimentally. The width can be accurately computed in some 
situations |571 1581 1ST] , although its derivation is, in general, non-trivial. 



4.1 V + A correlator in resonance theory 

At large Nc, the spin-1 correlators are meromorphic functions characterised by the position and 
residues of their poles. The V + A spectral function is given by 

hmU LL (tf c ^ = F 2 5(t) + lim £ Ff 5(t - M?) , (34) 

/ I IVl '* OO r , 

00 M?<M2 
J — 00 



which generates the moment integrals given by the positive meromorphic functions 

A ^ = ^ + E (35) 



with Mqo denoting the mass of the highest resonance in the resonance theory, the constants Mj 
and Fj are the mass and coupling constant of the vector and axial- vector states at LO in 1/Nc and 
F is the pion decay constant. In the limit when — » oo, this upper mass acts like an ultraviolet 
cut-off of the large Nc infinite resonance summation [HOI ESI EHl EZ1 • Taking the Q 2 — ► oo limit is 
not trivial at all since there is always an infinite set of resonance above any considered Q 2 , whose 
effects are in general not so clearly negligible. From now on, the limit — > oo will be always 
implicitly assumed. 



4.1.1 Perturbative and non-perturbative contributions in pQCD 

In this section we make an analysis of pQCD which will be relevant for the resonance calculation. 
We study how the pQCD amplitude is recovered in the deep euclidean thanks to a spectral function 
that behaves like ^ImIl LL (t) P at high energies. 



We will consider an auxiliary spectral function with the form 



\mU LL (tC C t D ee hmU LL (tf QCD ■ 9(t-t p ), (36) 



where the low energy region has been removed. It yields the moments 

, (n) f/ ^2\PQ CD f°° dt Q 2n l TmTT U ^PQCD 

tp (t + Q 2 ) n+1 n " 



In the deep euclidean, these moment integrals recover the pQCD moments A^(Q 2 ) pQCD with the 
proper a s running up to subdominant 1/Q 2m power corrections. 
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Figure 5: V + A Adler function from the pQCD space-like calculation up to 0(a s ), where the 
band represents the O (a 2 ) uncertainties. It is compared with A LL (Q 2 ) P ^ D , derived through dis- 
persion relations from considering the spectral function ^Im\l LL (t) pQCD ■ 9(t — t p ) (dashed curve, 
t p = (1.25 GeV) 2 and t p = (1.45 GeV) 2 for the upper an lower curves respectively). The dashed- 
dotted curve includes also the contributions from the pion, p(770) and ai(1260) whereas the solid 
lines consider instead the experimental data 22 for the spectral function at t < t p . 



The corresponding Adler function up to 0(a s ) is plotted in Fig. ©. The solid band covers 
an uncer tainty of ±(^) 2 . The value of t p was varied between the values t p = (1.26 GeV) 2 and 

t p = (1.45 GeV) 2 . However, the asymptotic behaviour A LL {Q 2 ) pQ °° = A LL {Q 2 ) pQ ° D + <D(a 2 s ) 

is only achieved at energies Q 2 S> 100 GeV 2 , by far much larger than the expected Q 2 ~ 2 GeV 2 . 
These discrepancies are much easier to understand through the observance of the Adler function at 
O(o2): 



A l (Q 



2 ,pQCD 



N C 
12vr 2 



Q 



+ 0(a s ) 



N C 
12vr 2 



1 



Q 2 + 



+ 0{a.) 



(38) 



where the term generates deviations of the order of 1% at Q 2 
coming from the 0(a 2 ) corrections. 



Q 

100 GeV 2 , larger than those 



This points out that although the high energy logarithmic dependence of Imll^ (t) P p ^ t D provides 
the proper asymptotic behaviour of the pQCD moment integrals, the non-perturbative region owns 
crucial information about how fast this limit is reached. If one includes the contribution of the 
resonances laying on t < t p (Goldstones from the chiral symmetry breaking, p(770) and ai(1260)) 
the asymptotic behaviour is again reached within the expected range of energies. Alternatively, one 
may consider the experimental data for ImII Li (t) at t < t p [22], getting similar results due to the 
moderate size of the width of these states (MjTj < -jy— ■ M 2 ). For Fig. © we have use the inputs 
M p = 0.77 GeV, M ai = 1.26 GeV, F = 92.4 MeV, F p = 154 MeV and F ai = 123 MeV gUiS]. 
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4.1.2 Perturbative and non-perturbative contributions in QCD^ 

We analyse in this section the conditions to recover pQCD in the deep euclidean region of 
momenta through a resonance theory. pQCD will naturally arise when summing up the infinite 
series of resonances, being the higher dimension contributions in the OPE a mere remnant from the 
discrete summation of poles. 



First of all, it becomes necessary to split the resonance series into two sub-series, 
U(-Q 2 f c ^°° = An(-Q 2 )2f~+An(-Q 2 )^:~: a perturbative part (denoted as An^T") 



where the resonances already lay on the perturbative QCD regime, with M 2 > M 2 ~ 2 GeV' 



2 



which will provide the asymptotic behaviour of the moments; and a non-perturbative part (namely 

f 2 < M 2 V , 



AIT c ^°°) with the resonance masses laying within the non-perturbative regime of QCD, M 2 < M 2 



which will drive how fast the moments converge to the pQCD description. M p is the mass of the 
last multiplet included in AU Nc ^°° . We will assume that the masses M,- are on increasing order, 
Mi < M 2 < ... 

N p F 2 Q 2n 

The non-perturbative sub-series A A" (Q 2 ) c ^°° = ; , r0 3 — rr is finite and it can be 

^ LL ^ 'non-p. Z-^l ( j\/f • + Q ) ™+ 

3 1 3 

M 2 

analytically expanded in powers of —4-, with Mf < M 2 ~ 2 GeV 2 . Hence, the pQCD behaviour 

Q 

and the logarithmic a s corrections are only recovered through the perturbative part of the series, 
which, in addition, generates extra contributions to the O i^ ^lm ^ OPE terms. 

We need to transfer the information from ^ImIl LL (t) pQ ° D to the discrete summation of infinite 
terms in QCD^. At this point one needs to make an assumption on the asymptotic structure of 
the mass spectrum at high energies. We will consider some smooth M 2 = f(j) dependence for 
j > p + 1. The pQCD spectral function is discretized in the perturbative region [M 2 , oo) through 
the step-like function H(t): 

OO -i 

H{t) = £ [9(M 2 -t)-e(M 2 ^- 1)} - lmU LL (M 2 f CD , (39) 
j=p+i 

providing the step-like interpolation H{M 2 ) = ^ImTl LL (M 2 ) pQCD for any j > p + 1 and H{t) = 

for t < M 2 The step function is defined as 9(y) = < ? p° r ^ ^ ^ . 
' | 1 for y > 

The difference between the moments A^(Q 2 ) H of this function and the original ones 
•^ll ( ( 5 2 )p^ D IS Si ven by the expression 

-»2n 



dt Q n . . (n) tr ^2\pQ CD 



*zm = i M2p (t ; Q i r+1 m) = az(qx::-[i-k\ m 
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with A^ L (Q 2 ) H > and A^ a number in the interval 



< A„" < lft|a - ( 5»^ 



Q 



Q 



(n + 1) — f- + In ^ + ° 



Q 



M; 



. Q 4 



l + 0(a s (M 2 )) 



(41) 

for n > 1, being (5M 2 =max{(5M ? 2 }°^ p+1 . This A n vanishes for any n faster than any logarithmic 
pQCD correction at Q 2 — ► oo. 



The last step consists on converting the step-like function H(t) into a narrow-width spectrum 
through the prescription 



3 y 3 ' 3 

for any j >p+l, which produces the moments 



ri(M 2 ) ■ 5M 2 = hmU LL (M 2 ) pQCD ■ 5M 2 , 



1 - A 



n+1 



(42) 



(43) 



with A n c a number in the range 



0. 



5M 2 



Q 2 +M 2 +SM 2 



, vanishing faster than any log. 



Hence, for high energies M 2 < Q 2 < Mxn the contribution coming from the perturbative part 
of the series takes the form 

» / M LL . pert. 

00 A(0 ) 

^ LL (-Q 2 Z c r r = n LL (-Q 2 Y QCD + E nt! , ( 44 ) 

m=l t * 

where, in addition, to the pQCD result one gets j> m power terms which vanish faster than a s (Q 2 ). 
The extraction of A( O ) per ' cannot be done by trivially expanding the resonance propagators in 

powers of since there are always infinite states with masses Mj| > Q 2 . They can be analytically 
computed just in some cases. In this work, they are recovered through numerical simulations. 

On the other hand, the non-perturbative part of the series produces just 1/Q 2m power terms at 
M 2 <^Q 2 : 

- a(o ll r n ~ p ' 

AIT J-Q 2 f c ~°° = Y V {2m) ' , (45) 

m=l ^ 

with a( o L ; m) r n - p - = (~ir zu F ! ( M ! r -1 - 

Matching the sum of the perturbative and non-perturbative sub-series with the OPE yields the 
constraints 

<<?;;> = A^r- + f 2 + t^i = o, 

p J=1 (46) 

= A(o;;r- - e^ 2 ^I = Ti^ G ^r> > o. 
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A priori one may think of two kind of structures for the spectrum: one may have a chaotic 
spectrum where the values of the masses {M 2 }^ =1 do not follow any particular law; or we may have 
an ordered spectrum where the masses M 2 are ruled by some asymptotic expression when n — > oo. 
Quark model and Regge theory studies seems to point out most likely the second option. In this 
case, once an asymptotic structure of the spectrum M 2 = f(n) is assumed , the prescription for the 
coupling F 2 from Eq. (|42j) is shown to be enough to recover the pQCD result. If one assumes besides 
that the couplings F 2 also behave smoothly for high values of n, then this condition becomes also 
necessary up to power corrections in Eq. (fl2|) : 



F 2 = 5M 2 



" „ 



(47) 



Actually, a variation of the choice from ilmll^ (M 2 ) pQCD to ^lmIl LL (M^_ 1 ) pQCD in Eq. (|4*7|) shows 



that both forms are equivalent since they differ by a term O ( aj ^.f )■ If the logarithm running 



ri 



of ^ImII ii (t) pQCD in Eq. (|47|) is replaced by a different log, then the corresponding deep euclidean 
amplitude follows a completely different logarithmic behaviour to that from pQCD. 

It is also possible to study the vector and axial correlators as separated entities. One 
could consider uncorrelated spectrums, such that they follow different asymptotic behaviours. 
There would be two different expressions 5M V , and 5M\, for the squared mass inter-spacing 
between vectors and axial- vectors, respectively, for high values of the masses. This would imply 
that, instead Eq. (|4*7|) . the couplings would be given by Fy ~ $M Vn ■ ^Imll vv (M Vn ) pQCD and 
F\ n ~ 5M\ n ■ llmU AA (M 2 A J QCD , in order to recover U vv (-Q 2 ) pQCD Uli AA {-Q 2 f CD at high 
Q 2 . The current analysis corresponds to the case where the vectors and axial- vectors follow a similar 
law, although it can be extended to more general frameworks in a straight-forward way. 



4.2 V — A correlator 



The left-right correlator and their moments are given in QCD^ by the limit: 

-lmU LR (tf c ^ = -F 2 5(t) + lim £ F 2 [—itj] S(t — Mj) , (48) 

which provides the moments 

A ^(Q 2 ) C -% + J* E t^tWt' (49) 



J 



with iTj the parity of the j-th multiplet, (—1) for vectors and (+1) for axial- vectors. This expressions 
are provided by the sum-rule integration in the circuit C = Ci n + C ou t shown in Fig. 



A (Q 2 ) = lim i — o xi U r Jt) = lim / — - „ -Indict), 

oo 7 o (t + Q ) n+ vr 

(50) 
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Figure 6: Integration circuit for the moment integrals. 

where the integration in C ou t is assumed to vanish as — ► oo. 

The V — A Green function will be also split into a perturbative and non-perturbative sub-series, 
where the couplings and masses where formerly fixed in the V + A analysis. In this case, the 
summation is not positively defined and one needs to specify the parity iTj of every state. We will 
consider a spectrum of multiplets with alternating parity where the lowest state is a vector, this is, 
TTj = At high euclidean momenta M 2 <C Q 2 <C M^, the perturbative sub-series shows the 

structure 

a / ss^LR , pert. 

00 AiO ) 

An M (-Q 2 cr = e , ( 5i ) 

m=l *■» 

where no dimension-zero term arises since both the vector and axial correlators reproduce the pQCD 
expression. The non-perturbative part of the series can be analytically expanded without further 
problems and the summation of both contributions must match the OPE through the constraints 

Mo;«r rt ' — f 2 + ti-^ F l = o, 

A<<;r' - ti-^FfMj = 0. (52) 

3=1 

3=1 

4.3 Study of some hadronical models for the spectrum 

In this section, we will study some of the current models for the meson spectrum. The res- 
onances with masses in the low energy region {M 2 < 2 GeV 2 ) are susceptible to stronger devi- 
ation with respect to the expected asymptotic behaviour M 2 = f(n) for high values of n. Be- 
low 2 GeV 2 , only two resonance multiplets are found, corresponding to the vector p(770) and 
the axial-vector ai(1260) |65j . The mass spectrum M 2 with n > 3 is interpolated through the 
first two multiplets above 2 GeV 2 , namely p(1450) with M 2 = M 2 , ~ 1.45 2 and ai(1640) with 



(O 



(2) 



(o LR ) 

X (4) ' 



(O 



(<■>) 
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M 2 = M 2 , ~ 1.64 GeV 2 jHSJ. We will include in the perturbative sub-series the multiplets with 
n > p = 2, leaving M 2 = M 2 and Mf = M 2 X for the non-perturbative part of the resonance series. 

In order to obtain the contribution to the condensates from the perturbative sub-series we 
will analyse the V + A and V — A Adler functions from All ^ °° minus the pQCD amplitude, 

(AA(Q 2 )^°° - A{Q 2 ) PQCD ). We will work only up to 0(a s ) in the pQCD distribution so 0{a 2 s ) 
uncertainties are expected in the derivation of the condensates. Considering the Adler function 
avoids complications like renormalization scale dependences or absolute convergence of the resonance 
series. Nevertheless, the V — A measurement is utterly improved by a cross-check analysis of 

n M (-Q 2 ). 



8 m AlO ' pert ' 



We perform the interpolation (AA(Qj)Z7 - A{Qf) ) = 2J — through 

m = l Ql 

eight points Qj = Q\ + AQ 2 (£ - 1), with Q\ = 4 GeV 2 , AQ 2 = 0.025 GeV 2 and I = 1,2.. .8 
From these eight constraints we extract the perturbative contribution to the first eight condensates 
A(0 (2m) . The interpolating point are varied to Q\ = 2, 3, 5, 6 GeV 2 in order to check the 
consistence of the procedure and the size of the uncertainties. It is not possible to go to higher 
or lower energies since, respectively, on one hand, the 0(a 2 ) corrections become comparable to 
the 0(1/Q G ) terms and, on the other, the OPE series breaks down within the very low energy 
range. The central values of the results represent the output for Q 2 = 4 GeV 2 whereas the error 
bars express the spreading coming from the analysis at different energies, providing an estimate 
of the 0{a s ) corrections. The errors in the condensates of dimension two, four and six coming 
from neglecting the 1/Q 2m terms with dimension beyond 2m = 16 and from the duality violating 
terms are expected to be rather suppressed. We will make some considerations about the duality 
violations in the analysis of the five dimensional model although a more exhaustive analysis can be 
found in Ref. HQ. 



4.3.1 QCD^ in 1 + 1 dimensions 



In the pioneer work on large Nc |12j . QCD was study at Nc — * oo in a configuration space of 
dimension 1 + 1. In this case the spectrum could be solved and followed an asymptotic behaviour 
with constant squared mass inter-spacing, SM 2 = 5A 2 . A similar type of spectrum appears as 
well in different QCD models [SOI EH E3- We will consider M 2 = A 2 + n5A 2 for n > 3, with 
5 A 2 = (Mf, - M 2 ) and A 2 = M 2 - 2,5 A 2 . 

If one remains at O(a s ) the imaginary part of the pQCD correlator is just a constant and 
therefore all the couplings are equal to F 2 = ■ 5A 2 ~ (122 MeV) 2 . Thence, the perturbative 
sub-series contributions to the condensates derived from the numerical analysis are 

A(0 ( "r rt - = (-45.8161°$?) • 10- 3 GeV 2 ~ -5F 2 , 
A(0 ( "r rt - = (41.00l°; 22 ) • 10- 3 GeV 4 ~ 8F 2 M 2 , 
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A(CV 
A(0 

* 

A(0 



LJ? > pert. 
2) ' 

Li? > pert. 
(4) ' 

Li? > pert. 
(6) ' 



(7.4349±8:g8S5) • 10- 3 GeV 2 
(-13.43lg;S|) • 1(T 3 GeV 4 
23.41?;?) • 10- 3 GeV 6 



-3F 2 M 2 , 
8 F 2 Mp . 



(53) 



In this model, one may actually compute the exact value of the condensates coming from the 
perturbative sub-series [2~TT 153] : 



A(0 { 



LL . pert . 
2m) ' 



-1 



A(0 



Lit . pert . 
(2m) ' 



m 
-I] 



12vr 2 

N c 



(5A 2 ) m B, 

2 f 2 ^) 

m 247r 2 v > 



in \SA 2 
B 



M 2 



2 5 A 2 



B„ 



2 5 A 2 



(54) 



being 

B 3 (x) 



B m (x) the 

ry ^ 3, /y>2 _|_ _1 ™, 



Bernoulli polynomials (-Bi(x) = x - 
..). For our choice of parameters one finds 



l 

2' 



Bo(x) 



X + 



6' 



iii , pert. 



12vr 2 



M 2 



1 



3A 2 
p' 



A(0 



(4) 



1 JVc 

2 12vr 



? m;, 



A(0 (2) ) 



LJ? . pert. 



N c 
24vr 2 



5A 2 , (5A 2 ) 2 

IT 2 ? 2M 4 , 
p p 



5A 2 



A(0 



Nc 



(4) 



A(<;r- 



24vr 



7 <5A 2 M 



<5A 2 
2M 2 



<5A 2 M 4 



24vr z 



<5A 2 

m 2 : 



-45.821 • 10~ 3 GeV 2 , 
41.08 • 1(T 3 GeV 4 , 

7.4357 • 1(T 3 GeV 2 , 
-13.45 • 1(T 3 GeV 4 , 
23.7 • 10~ 3 GeV 6 , 



(55) 



in total agreement with our former numerical calculation. This provides a check of the accuracy 
that will make our next numerical calculations reliable. The V — A condensates depend linearly 
on 5M 2 = 5 A 2 and both A( O ) per and A( O ) pcr have a strong dependence on the mass of 

J > (2m) ' > (2m) ' ° 1 

8M 2 

the lowest state included in the perturbative sub-series, M i. In the limit -j-pf — > 0, one recovers 



exactly the pQCD expression n(-Q 2 )^ D at 0(a°) from Eq. ® , with A( } pert ' ■- 

and A(0 L ^ ) pe " = 0, being t p = M 2 ,. The condensates result naturally of the expected size 
without imposing further constraints, pointing out the close relation between the different QCD 
scales, hadronic masses, mass inter-spacings, couplings and condensates. 



-1) 



m p 



We will take F = 92 .4 MeV, M p = 0.77 GeV, the asymptotic spectrum M„ = A 2 + nSA 2 for 
n > 3 and the pQCD correlator ^ImIl LL (q 2 ) as inputs, and the parameters F p , F ai and M ai will be 
left as free parameters. They can be recovered by demanding that the condensates ( O f 1 ), ( ) 
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Table 1: Contributions from the perturbative sub-series to the condensates within different hadron- 
ical models. The columns with Converging WSR1 and WSR2 corresponds to squared mass inter- 
spacings of the form 5M% ~ 1/V™ an d ~ 1/n, respectively. All the amplitudes are considered 
up to 0(a s ). 

and ( ) are zero. For a general resonance model, Eqs. ()46|) and ()52() yield the relations 
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In our case at 0(a s ) one finds 
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(56) 
(57) 



p ~ F ai ~ Fj> 3 , falling both res- 



The couplings remain of the order of the asymptotic value F, 
onance couplings and axial-vector mass within the expected range of values obtained in former 
phenomenological analysis |41U57llBlfll65| , These light states, laying by the non-perturbative QCD 
regime, suffer slight deviations from the asymptotic behaviour in such a way that the OPE is exactly 
recovered. 



The short distance matching fixes all the parameters in our analysis, providing a prediction for 
the condensates of higher dimension through Eqs. (|46|) and Q52|) : 
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They already fall pretty near the usual determinations ( O 
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The calculation can be taken up to 0(a s ) so the a s (Q 2 ) running in A^ L {Q 2 ) is recovered. 
The next-to-leading order computation is relevant in order to improve the determinations of the 
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Table 2: Predictions for the two first resonance multiplets -p(770) and ai(1260)— within each model 
for the inputs F = 92.4 MeV and M p = 0.77 GeV. We consider the amplitudes up to 0(a s ). 

condensates ( O ) and ( O ) , checking the impact coming from the perturbative QCD corrections 
in a s . The contributions to the condensates at this order that come from the perturbative sub-series 
may be found in the first column of Table lfl|). where we have considered the same values for a s (Q 2 ) 
and the pQCD correlator as in Section 3. In Table (J2J), it is possible to find the corresponding p(770) 
and ai(1260) couplings and masses, and the value of the condensates ( {4) ) and ( O ) derived from 

the short distance matching ( ) = ( O ) = ( O ^ ) = 0. Small variations arise with respect 

to the 0(cP s ) result, as expected if the asymptotic behaviour of resonance parameters depends 
smoothly on a s . However, the 0{a 2 s ) uncertainties increase the error in the V + A condensates by 
two orders of magnitude, where ( ) becomes now negative, pointing out the necessity of working 
at that order if more accurate determinations are required. 

4.3.2 Five dimensional spectrum 

Another available scenario that has appeared recently is the one provided by models in five 
dimensions |31l 132 j . Here the resonances appear as Kaluza-Klein modes from the quantization of 
the momentum in the fifth dimension, producing a four dimensional effective spectrum with the 
dependence M 2 ~ n 2 , this is, 5M 2 ~ n ~ \/M 2 . We will use the interpolation M n = A + n<5A, with 
5 A = (M a i^ —M p i) and A = M p i — 35A. The corresponding contributions to condensates coming from 
the perturbative sub-series are shown in Table (0). Taking these values and the inputs F and M p 
one derives the resonance parameters and the predictions for the condensates shown in Table ©. 
One finds an acceptable value for the ai(1260) mass although the values of the resonance couplings 
go high above the usual determinations. The predictions for the condensates ( O ) and ( ) 

appear slightly shifted with respect to 1+1 QCDoq though the positive sign for (0 L ) is properly 
restored. 

It is important to remark that this is not exactly a five dimensional calculation but an effective 
study of its spectrum. The series of infinite resonances in the 5D-theory are not regulated in the 
way of our cut-off, yielding some extra features as the existence of infinite Weinberg sum rules [22] . 
In the strict five dimensional calculation one must handle the full Kaluza-Klein pile (including also 

the first two resonances) with its precise dependence. For instance, one has M 2 = ^ (n + 
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and F 2 = ^ • ^ (n + §) for the RSI metric in Ref. 021, being h 
the infrared brane. This exact structure of the spectrum generates a V 
condensates, 



n M (-Q 2 ) 



12vr 



cxp 



3vrQ 

2m: 



37r/(4M p ) the position of 
A amplitudes without 

(59) 



For the perturbative QCD range of euclidean momenta the exponential factor produce a huge 
suppression (0(1(T 4 ) already for Q 2 = 2 GeV 2 and 0(1(T 9 ) for Q 2 = 10 GeV 2 ). Variations on 
the lightest Kaluza-Klein modes due to 0{a s ) corrections (present analysis) or perturbations on 
the metric |32| result into the appearance of 1/Q 2m power terms. The duality violating term in 
Eq. (J59|) affects our numerical OPE interpolation beyond the dimension-sixteen condensate, so the 
determinations are still safe. Nevertheless, the exponential may become enhanced with respect to 
the 1/Q 2m terms when considering moment integrals A^ k \Q 2 ) or components B^ k \Q 2 ) for large 
values of k. Eventual contributions might produce observable modifications to the usual OPE 
calculations |21j . 



4.3.3 Converging Weinberg sum rules 



For a general spectrum and within our considered cut-off regularization, large Nq 

f°° 1 

QCD do not obey in general the two Weinberg sum rules (WSR), / dt — ImU LR (t) = and 

Jo TT 

f°° 1 

/ dt i — Imll (i) = 0, since the infinite summations of resonances are not convergent. The ful- 
J o n 

filling of these two WSR immediately implies a dominant behaviour U LR (— Q 2 ) ~ in the deep 
euclidean. 

One may consider the covergence of the two WSR as an available hypothesis for the model 
building. In this case, the summations 

/OO 1 00 

dt t m -lmU LR (t) = - F 2 8 m '° + £ [ 7Tj] F 2 (M 2 ) m , (60) 
TT J— J 

must be zero (and therefore converging) for m = and m = 1 |16j . This demands that 
F 2 n ^ and F 2 M 2 n ^ 0. If we desire just two Weinberg sum rules, then F 2 (M 2 ) m 
must not vanish when n — » oo for m>2. If the asymptotic value of the couplings is fixed through 
F 2 ~ 5M 2 • ±ImIl LL (M 2 ) pQCD ~ 5M 2 , one finds in Eq. that F 2 (M 2 ) m ~ 5M 2 • (M 2 ) m must 
vanish when M 2 — » oo just for m = and m = 1; the inter-spacing must tend to zero as n 
goes to infinity faster than 5M 2 ~ ~ and slower than <5M 2 e~ M ™. 

We will analyse both extreme cases through two different modelings of the spectrum. On 
one hand, we will consider the model WSR1, owning the spectrum M 2 = A 2 + ^fn 5 A 2 
with 5A 2 = (M 2 , - M 2 ,)/( v / 4 - \/3) and A 2 = M 2 , - a/3 5A 2 . This provides the limit case 

SM 2 ~ -4= ~ For the other limit case we consider the model WSR2, with spectrum 
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Ml = A 2 + 5A 2 In (n), with 5A 2 = (M 2 , - M 2 ,)/(ln (4) - In (3)) and A 2 = M 2 , - In (3) 5A 2 . This 
model owns the high mass behaviour SM 2 ~ ^ ~ e _M ™. 

The results for both kinds of spectrums up to 0(a s ) can be found in the third and fourth columns 
of Tables (0) and ©. The axial-vector mass goes beyond the usual range 1 GeV< M ai < 1.26 GeV, 
whereas the resonance couplings decrease up to unphysical values, what seems to rule out this kind 
of converging-WSR models. 



4.4 Truncated resonance theory and Minimal Hadronical Approximation 

The contribution to the OPE condensates coming from the high mass multiplets 
vary depending on the model for the hadronic spectrum. However, all of them re- 
produce pQCD and generate contributions to the condensates of the right size -the 
standard hadronical parameters F 2 and M 2 — without demanding any further fine tuning, 



LL Pert tR Pert ^ o(F 2 (M 

> (2m) ' ' > (2m) ' \ ' 



P 



The matching with the OPE can be improved through a more exhaustive scanning of the reso- 
nance parameters (e.g. A 2 and <5A 2 in the 1+1 QCD^ spectrum) and analysing possible corrections 
to the asymptotic dependence of M 2 . Likewise, the second vector multiplet, with M 2 , ~ 2 GeV 2 , 
lies by the border of the non-perturbative QCD region and it may suffer still sizable variations with 
respect to the asymptotic behaviour. A deeper study should consider the three resonances p(770), 
ai(1260) and /o(1450) apart, within the non-perturbative sub-series. 

Nonetheless, in many situations the inclusion of the full infinite pile of hadronic states is quite 
involved. The information about higher mass states is rather poor, so one cannot yield a precise 
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determination of the perturbative sub-series AI1(— Q ) ncrt = II (— Q ) + 

m=l 



pert. ~ LL \ t « I ^2^1 Q2m 



This forces to truncate the resonance tower and to consider just a finite number of states, those 
in Ail(— Q 2 ) n ^ n _^ j setting AII(— Q 2 ) p ^ rt °° = 0. The separation between perturbative and non- 
perturbative series is not clear and, formally, one may consider as many states as desired within 
AII(— Q 2 )^ n _^ ■ Moreover, in order to reproduce the right long/short-distance behaviour of QCD 
(X-PT /OPE) one needs to include at least a minimal number of multiplets in the resonance theory. 
This is denoted as Minimal Hadronical Approximation (MHA) [10]. Strictly speaking, MHA refers 
only to Green-functions which are order parameters although some ansate replace An(— Q 2 )^ rt °° 
by a pQCD continuum [38j l39j. 

In the MHA analysis of the V — A correlator, one needs to keep at least the Goldstones from 
the chiral symmetry breaking, together with the first vector and axial-vector multiplets: 

„ , „o-.N n ^oo F 2 F 2 F 2 / 5M 2 \ 

n LR (-Q 2 ) ° = ? + n P n 9 1 n + O . 61 

lrK ' Q 2 M 2 + Q 2 M 2 ai +Q 2 \Q J 

The terms O f^p-) from AIl LR (—Q 2 )^ rt °° are dropped out in MHA. Through a short-distance 
matching to the dominant power behaviour of the OPE (U LR (-Q 2 ) ~ l/Q 6 ), the values of the 
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lightest resonance parameters are constrained, getting the familiar relations [161 1421 H3*[ . 

F 2 P - Fl = F 2 , F 2 M 2 - F 2 ai M 2 ai = 0. (62) 

The couplings and masses are now substituted by some effective parameters (F p / ai — > F p i ai and 
-^p/ai — * M p i ai ) where the information of the perturbative sub-series is also encoded. They suffer 
corrections proportional to A( O h ^ y ert - w ith respect to their original values in the full theory and 
some information is lost since the number of hadronic parameters decreases considerably. For the 
experimental values M p = 0.77 GeV and F p ~ 154 MeV, coming from the decay p° — > e + e~ [41| . 

one finds the predictions F ai = ^Fj - F 2 ~ 123 MeV, M ai = F p M p /F ai ~ 964 MeV and 
F 2 ~ 

( ) = --4 r F 2 M i ~ -4.7 • 1(T 3 GeV 6 . 

\ (6) ' p2 P 

1 ~ 

A further analysis of the vector and axial form factors |42| leads to F p = y/2F ~ 131 MeV, 
F ai = F ~ 92.4 MeV, M ai =V2M p ~ 1089 MeV and ( & L * ) = -2F 2 M^ ~ -6 • 10" 3 GeV 6 . One 
must be aware of the intrinsic uncertainty laying in the truncated resonance theory (which must 
not be misleadingly taken as model dependence). Performing a short-distance matching for a wider 
set of matrix elements modifies the value of the parameters and one risks to reach inconsistencies 
between the different constraints, eventually requiring the introduction the next resonance multiplet 
in the MHA description. 

The difference between the QCD parameters in MHA and those derived through other 
techniques is just due to the absence of the perturbative part of the resonance series. The 
uncertainty in the MHA short-distance matching to OPE is given by the size of 
A(0^\ ) pert ' ~ O (6M 2 (M 2 l ) m - i y The smooth variation of the resonance parameters from MHA 
to MHA+// hints its close relation with the full large Nq resonance theory. 

The employment of this approximation has led to a very successful large Nc phenomenology [131 ■ 
study of vector, axial- vector and scalar form- factors [1211001 > determination of xPT couplings [^ 
HIES [13 [111 El, description of two-point Green functions E3 OH E21 SSI and three-point QCD 
Green functions [441 1451 H5j . In addition, once the theory is well founded at leading order in 1/Nq, 
unitarity imposes serious constraints on the next-to-leading order effects. Thus, the final state 
interaction admits a well defined description within the MHA framework [551 1571 1581 1591 1601 1611 
I62 | I63 | IM] . Likewise, some calculations have affronted the problem of the renormalization and the 
radiative corrections to the couplings at next-to- leading order in 1/Nq [4*9"[ l5*U[ IBT ] l5*2 l l5*3*[ l5* H 15*5*]. 



5 Space-like region and local duality: resonance theory vs. OPE 

In this section we show how the resonance description matches the OPE in the euclidean domain. 
The 1+1 QCDoo model and the 5D-spectrum are considered although the first one yields the closest 
results to phenomenology. One can see on right-hand-side of Fig. (J7J) the comparison between the 
V + A Adler functions coming from the OPE and the large Nc resonance theory. We have plotted 
Q 6 Il LR (— Q 2 ) on the left-hand-side of Fig. (JJJ) in order to show the short distance 1/Q 6 behaviour 



26 



— f" 1 — 1—1 — r ~\~ l — 1—1 — r~| — r~r~i — r~p 



> 



o 

I 

""Si 
J 

tl ' 
o 



4 6 8 

I I I I I I I I I I I I I I I I I I I 



5D spectrum 
OPE 

— 1+1 QCD^ spectrum 
-- MHA 



Q 2 (GeV 2 ) 



1.4 f.*. 



1.3 



5D spectrum 
pQCD 

— 1+1 QCD^ spectrum 




Q 2 (GeV 2 ) 



Figure 7: V — A and V + A Adler functions in RxT^ 00 - 1 up to 0(a s ) and within MHA. They are 
compared to the results from the OPE and pQCD at that order. 



in a more transparent way. For the OPE and pQCD, we have taken the values from Section 3, with 
V — A operators up to dimension twelve. The MHA expression for Ii LR (—Q 2 ) is plotted for the 
values F p /V2 = F ai = F and M ai = V2M p . 

Although the 0(a 2 ) uncertainties are still sizable and subleading 1/Nc effects need to be anal- 
ysed, the euclidean amplitudes are already sensitive to the different structures of the spectrum. The 
1+1 QCDoo spectrum recovers quite accurately the OPE and pQCD. The amplitude coming from 
the five dimensional model does not agree for the choice of parameters in this work, pointing out 
the necessity of a further tuning of M p i and M a r in order to match the OPE at short distances. 
Finally, one can see that, although MHA produces &V — A correlator of the right order of magnitude 
and 1/Q e dependence, it does not completely agree the OPE determination. A better agreement is 
found if the short distance constraints of the form-factors are relaxed and only the two WSR are 
kept The inclusion in the MHA of the next resonance multiplet, corresponding to the p(1450), 
improves the p + a\ calculation and increases the accuracy of the determination 

Another way to study the properties of the amplitude in the euclidean region is through its 
moments at a fixed energy. In Section 3, we showed how the components B (Q ) suppress the 
dependence of the dispersive integral on the spectral function ^Imll(t) around t ~ Q 2 . The B k (Q 2 ) 
of the physical amplitudes oscillate as k grows and eventually damp off. The OPE is able to follow 
the physical B (k ' ' (Q 2 ) for the first values of k but, for a fixed Q 2 , it breaks down above some 
k OPE _ Hm . Increasing the energy or considering higher dimension operators allows increasing the 
range of validity of the OPE, k OPB _ lim . 

The large Nc resonance framework reaches a much further range, reproducing at the same 
time the good OPE results for k < k OPB _ Um . In Figs. (jHJ) and (jHJ), we can see the corresponding 
components from the experimental ansatz in the V + A and V — A channels compared with the 
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resonance description from the 1+1 QCD^ spectrum and the 5D spectrum, being again the first 
one clearly favoured. The situation is manifest in the V + A correlator, where the oscillation in the 
5D-model is much wider and one even finds opposite signs for the components of the two resonances 
models. In the V — A case, it is also compared with the MHA expression, with F p /^/2 = F ai = F 
and M ai = \/2M„. In general, the B k \Q 2 ) in the large Nq resonance theory follow a natural 
oscillating behaviour. It does not vanish for high values of k since the states own zero width, 
whereas the physical components eventually damp off due to the non-zero widths of the hadronic 
states. Nonetheless, through a careful estimate of the first components and their uncertainties, one 
is already able to distinguish between the different models for the spectrum and their couplings. 



4 — ♦ Exp. Ansatz 




Figure 8: Components B^ R (Q 2 ) of the V - A correlator for Q 2 = 2 GeV 2 and Q 2 = 10 GeV 2 . The 
experimental ansatz is compared to the large Nc resonance theory determinations and MHA. 



6 Time-like region and averaged spectral functions 



The QCD interaction at Nc — > oo is so strong that distorts the smooth pQCD spectral functions 
into a series of narrow-width resonances. However, the higher order corrections in 1/Nc provide the 
hadronic states with non-zero widths, so the physical amplitudes become again smooth. However, 
it is interesting to extract information from the experiment already at leading order in 1/Nq, so we 
will consider the averaged spectral functions introduced before in Section 2.3. Instead of working 
with the spectral function ^lmn(t) at a given positive energy t, it is rather convenient to employ the 
function a z (x), with the mapping x = ffipfj • a z{ x ) is then averaged by the distribution £ (x), with 
dispersion (Ax)| o = The average ^Imll(z) (a = (a z ) fa is essentially provided by the spectral 
function a z (x) in the interval \x\ < (Ax)^ a , i.e., by ^lmn(t) in the interval \t — z\ <2z (Ax)^ a . The 
interesting peculiarity here is that the averaged correlator only depends on the first a + 1 moments. 
As we saw in the former section, the first components B (z) are mainly governed by the large Nc 
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Figure 9: Components B^ L (Q 2 ) of the V + A correlator for Q 2 = 2 GeV 2 and Q 2 = 10 GeV 2 . The 
experimental ansatz is compared to the resonance expression at large Nq- 

resonance contributions, being the higher orders in k ruled by the subleading effects in 1/Nc- 

One may compare the experimental and theoretical {&z) ta pondered by distributions £ a (x) 
with a small enough. This procedure, suppresses the influence of the data away of the center of 
the distribution, t = z. However, if z lays near a resonance peak and a is taken too large, then the 
distribution £ a (x) may result too narrow and one has to consider the physical non-zero width of the 
hadronic state. 

This procedure could be useful for the analysis of the spectral functions from r-decays, where 
the data reaches just t < 3 GeV 2 . At low enough energies (for z < 1.5 GeV 2 and taking £ a (x) with 
a ~ 8), the region where there are no experimental data remains on the tail of the distribution 
£,a(x), yielding suppressed corrections. The range of energies around the p resonance is specially 
interesting, where very accurate data [221 I2H1 021 E51 [2E| and an exhaustive theoretical work on the 
resonance parameters already exists 021 E3 023 E3 EH 123 E2j ■ 

On Fig. ()1U|) . we can see how Rxr (oo) with the 1+1 QCD^ spectrum matches pQCD for 
z > 5 GeV 2 , where the slight discrepancy below demands a deeper analysis of the contributions from 
the non-zero dimension operators in the OPE. Actually, the fine agreement with the experimental 
ansatz around the /o(770) peak (z <1 GeV 2 ) points out again that the values in Table. ((21) for the p 
and a\ parameters lie on the proper range. For the set of parameters considered in this work, The 
5D spectrum shows again a slower convergence and large discrepancies with the experimental ansatz 
at z < 1 GeV 2 . MHA is also shown in the V — A case with the couplings F p /y/2 = F ai = F and 
masses M ai = y/2M p . At z < 1 GeV 2 , it provides an acceptable approximation of the experimental 
ansatz averaged amplitude. 

For energies beyond 1.5 GeV 2 the tail of the distribution £ a (x) reaches also the resonances with 
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squared masses around 2 — 3 GeV 2 and both MHA and the experimental ansatz with just r data 
loose reliability. This calls out for the inclusion of higher energy e + e~ data. A more exhaustive 
error analysis is also needed in order to yield a precise and accurate determinations of the p(770) 
and ai(1260) parameters at leading order in 1/Nc- 




1 1 2 3 

z (GeV 2 ) z (GeV 2 ) 



Figure 10: Averaged spectral functions ^lmll(z) = (a z )> for a = 8. 



7 Conclusions 

This paper explores the large Nq description of QCD through a theory with an infinite number 
of hadronic states. We show how it is possible to recover the OPE up to order a s in a systematic 
way, obtaining the corresponding a s (Q 2 ) running in H LL (— Q 2 ). Producing the precise anomalous 
dimensions in the condensates is still a hard task and the study at 0(a 2 ) is relegated to next works. 

The present analysis is based on three foundations. First, it is assumed that, in the deep 
euclidean domain, the QCD amplitudes and the different moments can be fairly described through 
the OPE. Second, the amplitudes in the large Nq limit accept a description in terms of an infinite 
exchange of narrow resonances that embodies the OPE; the correlators are meromorphic functions 
determined by the positions and residues of an infinite set of real poles. Third, in order to handle 
the infinite tower of hadronic states one has to assume a given asymptotic structure M 2 = f(n) for 
the spectrum at high energies and a smooth behaviour on n for the resonance couplings. Although 
we still lack a definitive theoretical explanation of the meson masses, we can nevertheless test (and 
eventually discard) some of the models currently considered. 1+1 QCDoo (M 2 ~ n) and the 5D 
theories (M 2 ~ n 2 ) seem to own the most favoured spectrums by the phenomenology, being the 
first one in closer agreement. 

The structure of the mass spectrum imposes constraints on the resonance couplings if we want 
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to recover the pQCD expression for the V + A correlator: 

Fl = Ml ■ ^U LL {Mlf CD + o(j^)\. 

The lightest mass states cannot be fixed through this procedure since pQCD breaks down and 
the asymptotic dependence of M 2 may suffer large a s corrections. The perturbative sub-series 
containing the high mass resonances is fixed by pQCD and our knowledge on M 2 . As the pQCD 
y+j4 correlator is recovered, the V — A amplitude shows an OPE structure of the form TL LR (— Q 2 ) ~ 
2 m ■ A final matching is performed by fixing the values of the lightest resonance parameters in 

order to ensure that H LR (—Q 2 ) goes like 1/Q 6 and that the first condensate in H LL (— Q 2 ) has 
dimension four. Although the deviations of F p and F ai with respect to the asymptotic values of F n 
are not large, they are essential to match RxT^°°^ and the OPE at order a s . The 1+1 QCD^ mass 
spectrum M 2 >3 = A 2 + n5A 2 yields the couplings 

F p = 139 ± 3 MeV , F ai = 138 ± 3 MeV , M ai = US0±\l MeV . 

The prediction for ( O J* ) agrees standard OPE determinations but there are still large uncertainties 
in the V + A sector where a deeper 0{a 2 ) study needs to be done. Considering pQCD just at 0{oP s ) 
produces small modifications on the resonance parameters, what ensures the stability of the whole 
procedure under a s corrections. 

In order to isolate the peculiarities of the quark-gluon and resonance pictures we develop and 
adapt a set of sum rules that are specially sensitive to the resonance features. In RxT °° , the 
combinations of moments B^ (Q 2 ) = J2i=i (Q 2 ) provided by the Legendre polynomials 

expose a characteristic oscillation which is also experimentally observed. The different models for 
the spectrum produce clearly different patterns of oscillation whose study can be used to discern the 
proper hadronical spectrum of QCD. The 1+1 QCD^ model reproduces the experimental estimate 
of the B ( \Q 2 ) up to k ~ 10 even for very low energies. For the V — A amplitude, the simple MHA 
expression provides a description as good as the one from the full RxT , pointing out the fact 
that it encodes an important portion of the QCD information. 

Another relevant technique for the determination of the resonance parameters is the employment 
of averaging distributions peaked around some energy t ~ z. Former works on Gaussian sum rules 
illustrate the averaging procedure It allows to isolate the data around a given energy region 
in order to study the parameters of a resonance laying in that interval or to perform checks of 
duality violations. In the present paper we have considered the family of distribution £ (x) with 
a power-like dependence, which tend to the Dirac delta when a — > oo. It is specially useful for 
the p(770) region where plenty of accurate experiments exist jUJ 1231 [2H E3 HB HZ| and the lack 
of data at high energies has little influence. A very good agreement between the experiment and 
the 1+1 QCDoq model is found for the averag ed amplitudes ±ImII(z) at z < 1 GeV 2 . The MHA 
provides as well a reasonable approximation. 

Through the splitting of the resonance series into perturbative and non-perturbative sub-series, 
we have understood how the truncated resonance theories -and more exactly MHA- connect the full 
large Nq theory. When the contribution to the condensates coming from the perturbative sub-series 
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is neglected, one introduces an error in the light resonance parameters. The terminology Minimal 
Hadronical "Ansatz" is then misleading since what we perform is simply an approximation with 
a clearly defined uncertainty. For the different models considered in this work, the contributions 



A( } ) pe " to the condensates are of the order of usual hadronical parameters, F 2 (M 2 )" 1 ^ 1 . The 
MHA parameters turn out to be of the same order of magnitude as in the full large Nc resonance 
theory, although even factor two discrepancies may arise, e.g., the value of the coupling F ai or 
the condensate (C ). This kind of considerations should be taken seriously when estimating the 
uncertainties in the MHA determinations. The inclusion of the next multiplet, p(1450), can improve 
the approximation in a stable way so the values of the condensates do no change drastically from 
MHA to MHA+p' 0, finding a smooth convergence from MHA to RxT 



To end with, I would like to comment some of the parallel work-lines arisen at the study of 
R X T V ■ As it was commented before, the B (Q 2 ) sum rules may be also used to fix the values of the 
condensates. They enhance both the low and the high energy regions where we have, respectively, 
experimental data and accurate theoretical descriptions. Likewise, by averaging the correlators 
through the distributions Ca(x), one is able to remove the influence of the condensates up to the 
desired dimension, leaving the pQCD contribution unchanged up to 0(a s ). An analysis of the 
average ^lmll vv (z)^ a from the e + e~ data would lead to an alternative determination of a s and 
high dimension condensates. 



The study under the OPE perspective has shown that, as far as the correlator shows a 1/Q 2m 
power-like structure in the deep euclidean, all the different functions present a similar OPE-like 
structure: 



n(-Q 2 ) = E 



(O 



(2m) 



Q 



2m 



A (n) (Q' 2 



B (k \Q 2 
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rn 

E 



(n,2m) 



Q 



2m ' 



(fc,2m) 

-imn 



^lmH(z) = E - 



(2m) 



IlmH(zf = E- 



1 _€a 

-imn. 



(2m) 



Up to 0(a s ), the averaged spectral function ^ImU(z) pQCD ' ia is just equal to the pQCD spectral 
function for any distribution £a(x). Possible duality violating terms in the QCD amplitudes could 
be eventually analysed in a similar way. 



All these techniques and considerations can be applied to other Green functions (scalar corre- 
lators, ...) and observables (form- factors,...). The analysis of exclusive channels would allow the 
extraction of the parameters of the resonance lagrangian at leading-order in 1 /Nc relevant for the 
different QCD matrix elements. These considerations are of particular interest in the case of more 
complicated Green-functions, e.g. three-point Green-functions, where the logarithmic behaviours 
are also originated by the infinite exchange of resonances through the different external legs. 
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